Renewable  and  Sustainable  Energy  Reviews  38  (2014)  572-580 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Renewable  and  Sustainable  Energy  Reviews 

journal  homepage:  www.elsevier.com/locate/rser 


Transmission  management  for  congested  power  system:  A  review  of 
concepts,  technical  challenges  and  development  of  a  new  methodology 


CrossMark 


Muhammad  Bachtiar  Nappua’*,  Ardiaty  Ariefa,  Ramesh  C.  Bansal 

a  Department  of  Electrical  Engineering,  University  of  Hasanuddin,  Makassar,  SULSEL  90245,  Indonesia 
b  Department  of  Electrical,  Electronic  and  Computer  Engineering,  University  of  Pretoria,  Pretoria,  South  Africa 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  27  June  2012 
Received  in  revised  form 
12  May  2014 
Accepted  24  May  2014 
Available  online  12  July  2014 

Keywords: 

Transmission  congestion  and  management 

Network  losses 

Economic  dispatch 

Optimal  power  flow 

Shift  factor 

Electricity  market 

Locational  marginal  price  (LMP) 


Transmission  networks  have  some  constraints  that  should  be  addressed  in  order  to  ensure  sufficient 
control  to  maintain  the  security  level  of  a  power  system  while  maximising  market  efficiency.  The  most 
obvious  drawback  of  transmission  constraints  is  a  congestion  problem  that  becomes  an  obstacle  to 
perfect  competition  among  the  market  participants  since  it  can  influence  spot  market  pricing.  As  the 
power  flow  violates  transmission  constraints,  redispatching  generating  units  is  required  and  this  will 
cause  the  price  at  every  node  to  vary.  This  manuscript  presents  concepts,  technical  challenges  and 
methodology  for  investigating  an  alternative  solution  to  the  redispatch  mechanism  and  then  formulates 
LMP  scheme  using  an  optimisation  technique  that  may  well  control  congestion  as  the  main  issue.  The 
LMP  scheme  are  varied  and  improved  to  take  into  account  the  energy  price,  congestion  revenue,  cost  of 
losses,  as  well  as  the  transmission  usage  tariff  by  utilising  shift  factor-based  optimal  power  flow 
(SF-OPF),  which  is  derived  from  the  well-known  DC  optimal  power  flow  (DC-OPF)  model. 
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1.  Introduction 

In  deregulated  electricity  markets,  transmission  networks  hold 
an  essential  role  in  supporting  interaction  between  producers  and 
customers  because  it  should  provide  an  unbiased  environment  for 
all  participants  [1],  Furthermore,  bottlenecks  in  transmission  lines 
are  an  obstacle  for  perfect  competition  among  market  participants 
[2],  One  obvious  drawback  of  transmission  constraint  is  conges¬ 
tion.  The  system  is  said  to  be  congested  when  both  parties  agree  to 
produce  and  consume  a  particular  amount  of  electric  power  but 
fails  to  do  this  because  of  a  transmission  network  exceeding  its 
thermal  limits.  This  eventually  brings  further  impacts  to  exercise 
market  power  that  can  cause  price  volatility  beyond  the  marginal 
costs  [3],  Therefore,  an  efficient  congestion  management  becomes 
one  of  the  most  important  indicators  of  the  design  of  a  power 
market  to  display  the  network  capability  in  a  competitive  market. 

The  fundamental  issue  in  a  competitive  market  is  the  market 
clearing  price  (MCP)  mechanism.  As  the  power  flow  violates  transmis¬ 
sion  constraints,  redispatch  of  generating  output  is  required  and  this 
causes  different  prices  at  every  node.  This  phenomenon  is  defined  as 
locational  marginal  prices  (LMPs),  also  known  as  load  pocket  [4],  Based 
on  these  facts,  the  relationship  between  MCP  and  transmission 
management  has  a  strong  relationship,  which  needs  to  be  assessed 
in  order  to  obtain  an  efficient  and  transparent  price  to  satisfy  all 
market  participants. 

The  DC-OPF  has  been  used  widely  to  primarily  manage  congestion 
through  computing  LMPs,  due  to  its  speed  and  strength.  Flowever 
when  the  LMP  scheme  takes  line  losses  into  account  the  advanta¬ 
geous  features  of  the  DC-OPF  based  LMP  mode  are  diminished. 

Therefore,  this  paper  examines  an  advanced  solution  for  the  re¬ 
dispatch  mechanism,  which  not  only  improves  the  computation  of 
congestion  and  losses,  but  also  formulates  a  new  LMP  scheme  using 
an  optimisation  technique  that  considers  congestion  as  the  crucial 
problem.  The  LMP  schemes  are  adapted  and  improved  to  take  into 
account  the  cost  of  losses,  congestion  revenue,  energy  prices,  as 
well  as  an  embedded  cost,  called  the  transmission  usage  tariff 
(TUT).  The  objective  is  to  support  the  development  of  a  standard 
market  design  in  managing  transmission  systems  which  promotes 
economic  efficiency,  lowers  delivered  energy  costs,  maintains 
power  system  reliability  and  mitigates  exercising  market  power. 
Accordingly,  three  schemes  of  the  LMP  are  introduced,  namely  LMP- 
lossless,  LMP-loss  and  LMP-TUT.  The  LMP-lossless  model  has  two 
components,  namely,  energy  price  and  transmission  congestion 
revenue.  The  second  scheme  in  this  LMP  modeling,  LMP-loss, 
includes  three  components:  energy  price,  transmission  congestion 
revenue  and  transmission  losses  cost.  LMP-TUT  is  formulated  based 
on  the  LMP-loss  but  takes  into  account  a  tariff  for  transmission 
usage  as  well.  These  schemes  are  developed  to  evaluate  the 
performance  of  the  proposed  method  of  shift  factor-based  optimal 
power  flow  (SF-OPF),  which  is  derived  from  the  DC-OPF  model. 

This  manuscript  is  organised  as  follows.  Section  2  elaborates 
issues  surrounding  transmission  management  in  the  electricity 
industry.  Section  3  gives  an  overview  of  economic  dispatch  and 


optimal  power  flow.  Section  4  explores  the  congestion  and  supply/ 
demand  equilibrium  and  Section  5  presents  a  generic  LMP  scheme 
for  congestion  management.  Section  6  proposes  an  improved 
method  of  OPF  for  LMP.  The  basic  tasks  of  the  improved  method 
and  its  conceptual  flowchart  are  presented  in  Section  7.  Finally,  the 
conclusion  of  this  paper  is  presented  in  Section  8. 


2.  Issues  surrounding  transmission  management  in  electricity 
industry 

Within  an  open  access  environment  transmission  management 
holds  a  vital  role  in  supporting  transactions  between  producers 
and  customers.  Bottlenecks  in  the  line  transmission  for  example, 
will  be  an  obstacle  to  perfect  competition  among  the  market 
participants.  Hence  the  operation  and  planning  of  a  transmission 
network  system  should  be  planned  in  an  effective  manner  [5,6], 
Fig.  1  shows  some  issues  faced  due  to  transmission  management  in 
a  deregulated  environment. 

One  obvious  drawback  of  transmission  constraint  is  congestion 
problems.  Congestion  is  a  result  of  transmission  constraints  limit¬ 
ing  network  capacity  that  interferes  with  power  transfer  from  a  set 
of  power  transactions  [7],  Two  other  significant  issues  that  should 
also  be  addressed  in  transmission  management  are  transmission 
usage  tariff  and  transmission  losses  [8-10],  Transmission  usage 
tariff  is  defined  as  embedded  cost  in  [6],  while  [11]  classified  it  as 
the  use  of  transmission  system  charge.  This  is  to  convert  standard 
operating  and  maintenance  costs  into  a  transmission  charge  cost, 
which  refers  to  the  previous  capital  cost  acquired  in  the  transmis¬ 
sion  infrastructure  development  and  maintenance.  The  last  aspect 
in  transmission  management  is  the  cost  of  losses.  Transmission 
losses  are  simply  defined  as  the  difference  between  the  total  power 
supply  from  generation  and  the  total  power  accepted  (demand)  by 
customers  in  the  system.  Even  though  the  impact  of  losses  may  be 
small  compared  to  other  potential  sources  of  market  inefficiency 
they  must  be  considered  as  well  [12],  To  meet  the  required  demand 


Fig.  1.  Issues  surrounding  transmission  management  in  a  deregulated  environment. 
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market  operators  have  to  make  sure  that  there  is  a  specific  supply 
for  network  losses.  In  terms  of  cost  minimisation  this  is  a  challen¬ 
ging  problem  to  find  an  appropriate  generator  responsible  for  the 
supply  of  losses  while  still  keeping  the  price  to  a  minimum.  Besides, 
in  the  planning  of  integration  of  renewable  distributed  generation 
(DG)  into  the  grid,  it  is  essential  evaluating  a  proper  assessment  on 
optimal  allocation  of  DCs  penetration  in  order  to  achieve  a  high 
efficient  electricity  market  through  minimising  total  energy  losses 
and  maximising  social  welfare  [13-15], 

Transmission  management  is  faced  dealing  with  many  problems, 
however  the  transmission  congestion  issue  remains  the  central 
problem  in  the  new  electricity  structure  [16,17].  The  objective  of 
managing  congestion  is  the  cost  of  re-dispatch  of  generation  as 
constraints  are  violated.  It  also  guarantees  that  there  is  a  sufficient 
revenue  to  cover  the  cost  of  transmission  system  operators.  The 
revenue  can  be  used  in  various  ways  to  address  congestion,  such  as 
to  promote  sufficient  transmission  construction  planning  which  is 
not  readily  done,  like  installing  new  generation,  expansion  of  net¬ 
works  or  employing  FACTS  devices  [18],  A  study  in  [19]  shows  that  by 
using  FACTS  devices  the  transmission  congestion  is  mitigated  with¬ 
out  disturbing  economic  matters. 

Independent  system  operators  (ISOs)  usually  observe  the 
transactions  and  control  the  state  of  the  system,  taking  part  in 
handling  the  network  congestion  management  [20,21].  ISOs  are 
being  challenged  to  develop  a  set  of  regulations  to  control  the 
security  level  of  power  systems  and  ensure  that  they  are  at 
acceptable  level  while  keeping  the  efficiency  of  the  power  market 
high  [22],  This  implies  that  market  operators  should  alleviate 
network  congestion;  maintain  the  security  and  efficiency  of  power 
system  operation  [23],  in  order  to  ensure  that  all  market  partici¬ 
pants  have  the  same  rights  to  access  a  transmission  system 
without  any  discrimination  [24].  Congestion  levels  typically  deter¬ 
mine  the  security  of  a  power  system  which  would  have  further 
consequences  on  market  transaction  and  energy  prices. 

3.  Economic  dispatch  and  optimal  power  flow 

3.1.  Economic  dispatch 

The  definition  of  economic  dispatch  as  cited  in  [25]  is  “the 
operation  of  generation  facilities  to  produce  energy  at  the  lowest 
cost  to  reliably  serve  customers,  recognising  any  operational  limits 
of  generation  and  transmission  facilities”.  While  reference  [26] 
defines  economic  dispatch  as  “the  process  of  allocating  generation 
levels  to  the  generating  unit  in  the  mix,  so  that  the  system  load 
may  be  supplied  entirely  and  most  economically”.  The  production 
cost  of  generation  is  analysed  during  the  dispatch,  subject  to 
data  which  is  concerning  fuel  cost  and  electrical  power  output. 
A  quadratic  equation  is  used  to  approximate  the  cost  function 
along  with  several  cost  coefficients.  The  key  objective  of  the 
economic  dispatch  problem  is  to  find  a  set  of  active  power 
delivered  by  the  committed  generators  to  satisfy  at  any  time  the 
required  demand  subject  to  the  unit  technical  limits  and  at  the 
lowest  production  cost.  For  this  reason  it  is  of  great  importance  to 
solve  this  problem  as  fast  and  precisely  as  possible. 

3.2.  Optimal  power  flow 

The  development  of  electricity  markets  all  over  the  world  has 
brought  new  challenges  to  the  scientific  community,  to  the  players 
and  to  the  regulatory  authorities.  One  of  these  challenges  is  the 
uncertainty  which  became  a  structural  element  in  new  environ¬ 
ment  [27].  All  market  participants  have  to  be  able  to  deal  with  it  to 
guarantee  the  appropriate  power  system  planning  and  operation 
as  well  as  its  own  economical  liquidity  [28].  A  well  accepted  tool 


used  for  controlling  and  assessing  the  appropriate  and  secure 
operation  of  a  power  system  is  optimal  power  flow  (OPF).  This 
economic  dispatch  technique  is  already  applied  in  many  energy 
management  systems  (EMS)  but  has  some  security  constraints. 
The  concept  of  OPF  was  defined  by  Carpentier  in  1962  [29]  as  an 
extension  of  traditional  economic  dispatch  of  power  to  resolve  the 
optimal  settings  for  control  variables  while  considering  various 
constraints  [30], 

The  OPF  comprises  mainly  of  an  economic  dispatch  mechan¬ 
ism,  a  power  flow  program  and  a  redispatch  generation  algorithm 
as  violation  of  security  constraints  occurs.  Such  security  con¬ 
straints  are  primarily  due  to  the  thermal  limits  of  transmission 
lines  and  capacity  limits  of  generation  [31],  The  economic  dispatch 
yields  a  single  incremental  cost  known  as  lambda  when  the 
security  constraints  are  not  violated.  Once  the  violation  occurs, 
redispatch  is  needed  and  economic  dispatch  will  produce  a  set  of 
different  lambda  at  different  nodes.  Lambda  generated  at  this 
stage  could  be  considered  as  a  marginal  cost  for  an  appropriate 
generating  bus  to  reschedule  power  supply  such  as  to  meet  all 
required  demand  of  the  bus  at  minimum  cost.  Marginal  cost  is 
referred  to  as  the  incremental  cost  subject  to  the  cost  of  producing 
additional  output  or  one  more  unit  of  output,  1  MW  h. 

3.3.  DC  optimal  power  flow 

The  optimal  power  flow  is  a  very  large  and  complex  mathema¬ 
tical  program.  An  OPF  can  be  described  as  the  minimisation  of  real 
power  generation  cost  in  an  interconnected  power  system  while 
real  and  reactive  power,  transformer  taps  and  phase-shift  angles  are 
controllable  and  various  inequality  constraints  are  required  [32],  Its 
procedure  consists  of  methods  of  employing  power  flow  techniques 
for  the  economic  dispatch  while  definite  controllable  variables  are 
adjusted  to  minimise  the  objective  functions  such  as  the  cost  of 
active  power  generation  or  the  power  losses  while  satisfying 
physical  and  operating  limits  on  various  controls,  dependent 
variables  and  function  of  variables. 

While  some  authors  [33-39]  have  used  the  AC  power  flow 
model,  others  [39-46]  have  used  the  DC  optimal  power  flow  model. 
The  AC  optimal  load  flow  problem  such  as  the  OPF  based  on 
Gradient  and  Newton’s  methods  consists  of  finding  the  active  and 
reactive  power  outputs  and  the  voltage  magnitudes  of  any 
generator  unit  in  order  to  minimise  the  operating  cost  while 
meeting  various  security  constraints.  Whereas  the  DC  OPF  approach 
is  applied  to  compute  the  load  flow  in  the  model  characterised  by 
ignoring  losses  and  the  focus  is  only  on  real  power.  The  DC  Optimal 
Power  Flow  (DC-OPF)  is  a  Quadratic  Separable  Program  with  the 
aim  of  resolving  the  interests  of  both,  suppliers  (who  are  wishing  to 
sell  energy  as  expensive  as  possible,  and  buyers  (who  are  wishing  to 
buy  energy  as  cheap  as  possible)  [47],  This  solution  is  widely  used  as 
this  is  the  conversion  of  the  ac  approximation  to  a  more  simple 
linear  circuit  analysis  problem.  Even  though  using  a  full  AC  power 
flow  is  the  most  precise  computation  [48],  however  owing  to  its 
complexity  and  non-linearity  features  it  may  obscure  parameters 
correlation.  On  the  other  hand  the  DC  OPF  methods  and  software 
are  simple,  its  models  can  be  optimised  efficiently  as  well  as 
minimal  network  data  required  is  fairly  easy  to  acquire  [42].  Thus 
the  DC  power  flow  approximation  is  preferred  and  has  been  used 
extensively  in  many  literatures  when  calculating  the  nodal  prices 
and  analysing  bottlenecks  over  transmission  lines  [40,49]. 

In  [29]  the  DC  power  flow  model  is  developed  with  some 
considerations  such  as;  voltage  magnitudes  are  assumed  to  be 
constant  at  every  bus,  the  different  phase  angles  of  bus  voltages 
between  two  buses  are  small  and  transmission  resistance  is  very 
small  compared  to  its  admittance,  therefore  there  are  no  losses  in 
this  model.  These  considerations  construct  a  model  that  is  a 
reasonable  first  approximation  for  a  real  power  system  which  is 
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only  slightly  nonlinear  in  normal  steady  state  operation.  The 
model  has  the  benefit  of  speed  of  computation  and  other  useful 
properties  such  as  [50];  linearity  feature,  this  means  the  flows  that 
are  directly  attributable  to  any  transaction  are  doubled  if  the  MW 
in  that  transaction  from  one  zone  to  another  is  also  doubled  and 
the  superposition  feature,  means  the  flows  on  the  networks  can  be 
divided  into  a  sum  of  components  each  directly  attributable  to  a 
transaction  on  the  system. 


4.  Congestion  and  supply-demand  equilibrium 

Social  welfare  holds  a  significant  role  when  assessing  the 
condition  of  market  economics  in  a  power  system,  especially  at 
the  supply-demand  equilibrium  under  congestion  circumstances 
[50,51],  Social  welfare  is  a  function  of  the  combination  of  the  cost 
of  energy  production  and  the  cost  of  customer  benefit  based  on 
their  ability  to  pay  for  it.  In  other  words  it  is  generally  a  function  of 
all  utilities  in  society  [52],  Both  supplier  g  and  customer  d  will 
have  supplier  surplus  SS(g)  and  customer  surplus  CS(d)  respec¬ 
tively  [53]. 

Nevertheless,  these  well-organised  systems  most  likely  become 
imperfect  when  congestion  occurs  over  a  certain  transmission  line. 
It  may  lead  to  market  inefficiency  due  to  the  thermal  limit  of  the 
network.  An  example  of  this  phenomenon  is  shown  in  Fig.  2  for  a 
two  zone  system  connected  by  a  tie  line  and  ignoring  the  line  loss. 
Let  both  zone  1  and  zone  2  have  a  constant  load  of  125  MW  where 
their  generating  capacities  are  alike  and  as  much  as  300  MW. 
It  is  assumed  that  all  generators  use  their  marginal  cost  for  their 
bidding. 

As  shown  in  Fig.  2(a),  in  order  to  maximise  their  own  profits 
each  supplier  involves  bidding  its  marginal  cost  when  a  generator 
is  a  price  taker.  Whilst  there  is  no  congestion  all  250  MW  of 
demand  will  be  directly  bought  from  zone  2  since  its  incremental 
cost  is  cheaper  than  zone  1.  The  total  purchase  cost  at  $  10/MW  h 
will  be  $2500/h.  Yet  if  power  flow  over  the  tie  line  is  limited  to 
75  MW  as  in  Fig.  2  (b)  only  200  MW  will  then  be  bought  from 
generator  2  at  $  10/MW  h.  Now  the  total  purchase  cost  will  be 
$3000/h  as  the  remaining  50  MW  will  be  ordered  from  generator 
1  at  $  20/MW  h.  It  is  obvious  that  congestion  has  created  an 
inefficiency  in  the  power  market  by  $500/h  or  equal  to  20%  of  the 
optimum  costs. 

In  the  above  case  no  strategic  bidding  from  the  generators  is 
employed.  Indeed  market  inefficiency  due  to  the  congestion  problem 
would  suffer  further  from  applying  strategic  behaviour  of  generators 
[54].  This  strategy  is  defined  as  an  effort  to  maximise  profits  due  to 
imperfections  in  the  market,  especially  when  congestion  exists. 
Strategic  bidding  is  generally  characterised  by  generators'  actions, 
which  bid  energy  price  not  at  their  actual  incremental  cost.  Once  a 

a 

300  MW  300  MW 


125  MW  125  MW 

Fig.  2.  A  simple  transaction  between  two  zone  systems  (a)  without  congestion 
(b)  congestion  with  thermal  limit. 


generator  can  advantageously  raise  its  profits  by  strategic  bidding  as 
a  logical  consequence,  market  power  will  be  created  [55], 


5.  Generic  LMP  scheme  for  congestion  management 

5.3.  Definition  of  locational  marginal  prices  (LMP) 

LMP  is  the  key  factor  to  identify  the  spot  price  and  to  manage 
transmission  congestion  [56],  The  LMP  is  utilised  to  calculate  power 
dispatch  schedules  by  maximising  the  social  welfare  function  [57], 
LMP  methodology  has  been  implemented  at  some  independent 
system  operators,  for  example:  California  ISO,  PJM,  New  York  ISO, 
ISO-New  England,  Midwest  ISO,  ERCOT,  etc  [58-60].  LMP  is  defined 
in  [53]  as:  “the  marginal  cost  of  supplying  the  next  increment  of 
electric  energy  at  a  specific  bus  considering  the  generation  marginal 
cost  and  the  physical  aspects  of  the  transmission  system".  In  other 
words,  the  LMP  is  the  cost  to  serve  one  additional  MW  of  load  at  a 
specific  location  using  the  lowest  production  cost  of  all  available 
generators  while  observing  all  transmission  constraints.  From  the 
viewpoint  of  generation  and  transmission  planning  it  is  fundamen¬ 
tal  to  always  calculate  and  forecast  the  value  of  LMP  which  may  be 
obtained  using  the  traditional  production  cost  optimisation  model 
[38,61],  However,  the  implementation  of  LMP  must  be  able  to 
manage  various  impracticalities  such  as  determining  prices  when 
systems  run  out  of  controls  and  transmission  constraints  are 
violated  or  there  is  an  excess  of  generation,  i.e.  total  generation  is 
greater  than  total  load  [59], 

5.2.  Approaches  of  DC  optimal  power  flow  model 

The  DC-OPF  model  is  a  common  optimisation  based  technique 
that  has  been  proposed  by  many  market  operators  and  is  being 
widely  used  in  various  ways  both  for  dispatching  power  and 
clearing  energy  markets  to  decide  the  LMP.  This  methodology 
has  become  the  leading  approach  in  electrical  power  markets.  The 
DC  optimal  power  flow  model  has  been  utilised  in  the  electricity 
industry  to  compute  LMPs  due  to  its  high  speed  of  convergence, 
simplicity  and  robustness  [41],  It  is  particularly  used  in  market 
simulation  and  planning  [62,63]. 

Generally  the  DC-OPF  is  used  for  security  constrained  economic 
dispatch  and  redispatch  when  controlling  transmission  congestion 
while  maximising  the  economic  power  transfer  capability  of  the 
transmission  system  without  violating  its  constraints  [64,65], 
These  auction  models  perform  system  security  by  means  of  simple 
power  flow  constraints  and  hence  can  be  stated  as  security  constraint 
OPF  problems. 

Ref.  [66]  gives  a  tutorial  review  of  the  use  of  optimal  power  flow 
approximation  to  calculate  LMPs  and  congestion  costs.  The  standard 
DC-OPF  scheme  for  LMP  calculation  can  be  modeled  as  the  mini¬ 
misation  of  the  total  production  cost  subject  to  energy  balance  and 
transmission  constraints.  This  may  be  written  as  follows: 

5.2.3.  Objective  function 

Let  the  total  production  cost  and  the  total  customer  benefit  are 
given  by: 

Ci(Pa)  =  ciiPa  +  bjPGi  +  q  ( 1 ) 

Bi(P  Di) =  djPpj+ ejPpj  + fj  (2) 

where  C,  is  production  cost  at  bus  i  ($/MW  h),  33,  is  customer 
benefit  at  bus  i  ($/MW  h),  Pa  is  active  power  output  of  ith 
generator  (MW),  PDi  is  active  power  demand  of  ith  load  (MW), 
a„  b„  cf  are  coefficients  for  power  production  cost  function  and  d„ 
ef,  f  are  coefficients  for  customer  benefit  cost  function. 
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The  objective  function  is  then  to  maximise  the  total  social 
welfare  (TSW)  which  also  equals  to  minimise  the  total  social  cost; 

max  Bi(PDi)~  £  C,(PCi)  (3) 

i  =  1  i  =  1 

where  nD  is  number  of  load  and  nG  is  number  of  generator. 

TSW  is  the  difference  between  total  customer  benefit  B,(PDi) 
and  total  production  cost  Cj(Pa),  which  is  going  to  be  maximised 
in  the  pool. 

5.2.2.  Constraints 

In  a  pool  based  dispatch  mechanism  an  optimal  solution  is 
required  subject  to  a  set  of  practical  constraints  as  follows: 

Nbus 

Pi  =  2  [M'S.-'Sj] 

i  =  1 
O'  #  0 

Pi  =  Pa~PDi  (j=\,...,Nbus) 

PT<Pa<Pcr  (i=  1.2,  ...Nc) 

Pi/S)  <  P|f x  (i=  1,2,  ...Nbus)  (4) 

where  P,  is  net  injection  at  bus  i,  Pjj" ax  is  upper  limit  of  power  flow 
between  busses  i-j,  by  is  susceptance  of  the  branch  connecting 
busses  i-j  (By  =  —(1/Xyj),  Sb  Sj  are  voltage  angle  of  bus  i  and  j  and 
Nbus  is  number  of  bus. 

5.2.3.  Optimality  condition 

From  the  objective  function  and  some  constraints  Lagrangian 
function  is  then  performed  after  selecting  a  reference  node.  An 
optimality  condition  will  result  in  a  set  of  decision  variables  and 
multipliers,  z,  which  may  be  written  as  follows. 

Z=[P  1  P2  ...Pn  S2...8„_i  A 1  A2  ■  ■  -An  P}  p2--Pm\T  (5) 

where  An,  pm  are  Lagrange  multipliers. 

If  there  is  no  transmission  congestion  it  means  there  is  no 
violation  on  transmission  constraints  and  no  redispatch  is  needed. 
In  this  case  the  economic  dispatch  program  creates  a  single 
incremental  cost  (lambda  or  Lagrange  multipliers),  modeling  the 
behavior  of  a  basic  economic  dispatch  algorithm.  Once  transmis¬ 
sion  congestion  occurs  security  limits  are  violated  and  redispatch 
is  required  to  meet  the  security  constraints;  the  result  is  a  set  of 
different  Lagrange  multipliers  at  different  buses  in  the  power  flow 
model.  The  optimal  price  is  represented  through  the  Lagrange 
multipliers  A.  The  Lagrange  multipliers  may  be  considered  to  be 
the  marginal  cost  or  incremental  cost  at  each  generating  bus  to 
redispatch  the  generation  in  a  manner  to  produce  energy  to  serve 
load  in  the  modeled  area  at  minimum  cost.  Congestion  cost  in 
these  models  are  formulated  as  ty=Aj—Ai, 

From  the  LMP  scheme  above,  it  can  be  seen  that  typically 
DC-OPF  is  developed  for  the  production  cost  model  given  the 
information  of  generation,  transmission  and  load  data.  The  natural 
fit  into  the  linear  programming  model  for  simplicity,  robustness 
and  efficiency  is  the  main  factor  why  a  DC-OPF  approach  is  the 
most  popular  one  to  be  used  in  performing  LMP  formulation. 
The  value  of  LMPs,  energy  price  and  congestion  cost  will  follow  a 
perfect  linear  model  with  zero  loss  prices,  in  particular  if  the  line 
loss  is  ignored  in  the  DC-OPF  model. 

Nevertheless  challenges  arise  if  loss  must  be  considered.  This  is 
due  to  the  non-linear  or  quadratic  relationship  between  line  loss  and 
line  current  or  line  flow.  These  weaknesses  are  quite  similar  even 
when  an  ac  optimal  power  flow  approach  is  employed.  Moreover, 
in  53]  it  is  recommended  to  include  the  cost  of  marginal  losses 
when  performing  the  LMP  model  as  it  is  known  that  resistance  in 
transmission  lines  is  inevitable. 

This  implies  the  LMP  model  should  encompass  generation 
marginal  cost,  congestion  cost,  as  well  as  cost  of  marginal  losses 


[67].  The  LMP  model  with  such  a  scheme  has  already  been  practiced 
in  many  ISOs  such  as  PJM,  CAISO,  NYISO,  etc.  In  [68]  the  authors  give  a 
solution  of  how  to  include  line  loss  when  using  the  DC-OPF  model  for 
economic  dispatch  mechanism.  This  solution  applied  to  [29]’s  work 
would  result  in  a  more  comprehensive  outcome. 

5.3.  Incorporating  line  loss  in  to  the  model 

The  line  loss  is  important  since  it  can  cause  major  impact  on 
optimal  economic  dispatch  depending  on  the  network  [69].  The 
losses  based  LMP  scheme  can  be  modeled  as  the  minimisation  of  the 
total  production  cost  subject  to  energy  balance  and  transmission 
constraints  [70],  Thus  LMP  incorporating  losses  comprises  of  three 
components:  marginal  energy  price,  marginal  congestion  price  and 
marginal  loss  price  [62,71,72],  Hence  the  first  constraint  in  Eq.  (4), 
bus  power  balance,  with  line  loss  taken  into  account  becomes: 

Nbus  Nf,us  rr;;-t 

Pi=  2  mst-Sji-  i  Um-Sj?  <6) 

i = 1  i= 1  LZ J 

0V0  0'#0 

Since  the  line  loss  in  this  formulation  is  a  function  of  voltage 
angle  differences  the  total  power  loss  does  not  depend  on  the 
selection  of  a  reference  node.  One  particular  voltage  angle  may  be 
chosen  to  become  zero  to  represent  a  reference  node  yet  the  overall 
angle  differences  will  always  remain  constant.  The  supply  for  power 
loss  is  distributed  among  all  generating  units.  While  distributing  the 
power  loss  the  angle  difference  over  congested  lines  must  be  kept 
same  to  satisfy  the  binding  constraint.  The  disadvantage  of  this 
technique  is  that  such  a  procedure  performs  a  uniform  amount  of 
total  losses  and  ends  up  in  a  high  cost.  Hence  an  alternative  solution 
is  needed  to  perform  an  optimal  LMP  model  which  discovers  a 
lower  power  loss  so  that  network  losses  cost  may  also  be  reduced. 

Therefore  in  the  following  section  a  different  approach  will  be 
introduced  by  using  the  shift  factor  method  with  the  intention  of 
maintaining  the  linearity  and  superposition  features  of  the  LMP 
model  while  this  will  be  able  to  optimise  for  both  congestion  cost 
and  losses  cost. 

6.  An  improved  method  of  OPF  for  locational  marginal  prices 

6.1.  Basic  concept  of  shift  factor 

From  a  power  flow  perspective  a  transaction  can  be  seen  as  an 
amount  of  power  that  is  injected  by  a  generator  into  the  system  at 
one  zone  and  received  by  a  load  at  another  zone.  A  DC  power 
flow  owns  a  linearity  characteristic  that  can  be  occupied  to 
calculate  the  transaction  amount  that  would  provide  rise  to  a 
specific  power  flow,  such  as  an  interface  limit.  The  coefficient  of 
the  linear  correlation  between  a  transaction  quantity  and  the  flow 
on  a  line  is  called  the  shift  factor  (SF).  The  shift  factor  is  the  actual 
component  that  determines  the  power  flow  over  a  given  transmis¬ 
sion  line  from  the  source  node  (generation)  to  the  sink  node  (load) 
[73],  It  is  characterised  by  four  attributes;  a  reference  node,  a 
particular  node,  a  particular  line  with  reference  direction  and  the 
value  of  the  shift  factor.  Shift  factor  is  also  known  as  power 
transfer  distribution  factor  (PTDF).  The  value  of  the  shift  factor  of 
line  /  with  respect  to  bus  i  is  related  to  the  amount  of  one  change 
of  active  power  flow  in  a  reference  direction  over  a  given 
transmission  line  /  to  another  change  both  in  withdrawal  at  the 
reference  node  and  in  injection  power  at  bus  i.  In  other  words,  the 
shift  factor  can  be  seen  as  the  fraction  of  transaction  amount  in 
line  l  due  to  an  injection  change  at  node  i. 

A  flow  in  line  / 

Sr  =  —. — : — 7 - : - : :  (/) 

A  injection  at  node  i 
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6.2.  Formulation  of  shift  factor 

Basically,  the  shift  factor  can  result  from  three  stages  based  on 
DC  load  flow  network  model,  as  follows  [73]: 

•  Stage  1:  Making  sensitivity  equation  of  phase  angles  as  the 
function  of  node  injections. 

•  Stage  2:  Making  sensitivity  equation  of  branch  flows  as  a 
function  of  phase  angles. 

•  Stage  3:  Substituting  the  result  of  stage  1  into  stage  2. 

7.  Basic  tasks  of  the  improved  method  and  conceptual 
flowchart  of  the  schemes 

The  methodology  of  this  research  comprises  of  several  basic 
objectives  and  is  divided  into  a  number  of  steps  as  follows. 


In  order  for  the  congestion  cost  and  locational  marginal  prices 
to  obtain  an  optimum  solution,  the  first  order  optimality  condi¬ 
tions  of  Karush-Kuhn-Tucker  (KKT)  should  be  met,  which  are 
derived  from  the  above  Lagrange  function.  This  results  in  loca¬ 
tional  marginal  prices  for  both  suppliers  and  customers: 


Pi  —  ^i  +  Z  Pmin.flow.m  ■  2  Pmax.flow.m  ■ 


(10) 


such  that 


dflowm  c  dflowm  c 
d(Pc)  ~  m’  d(PD )  ~  m 


7.3.2.  Scheme  II:  Named  as  LMP-loss  model  includes  three 
components:  energy  price,  transmission  congestion  revenue  and 
transmission  losses  cost  (LC) 


7.1.  Identification  of  primary  data 

The  identification  of  primary  data  is  divided  into  three  data 
collection  tasks.  The  first  task  is  generation  profile  data  including 
production  cost  parameters  to  perform  a  generator's  marginal 
bids,  location  of  generators  and  a  generator’s  output  limits.  This  is 
continued  in  the  second  task,  namely  identification  of  power 
system  network  profile  data  such  as:  number  of  branches,  resis¬ 
tance  and  reactance  values  of  transmission  lines,  thermal  limits  or 
maximum  capacity  of  each  branch  and  network  charge  rate  of  the 
branch.  Identification  of  demand  profile  data  is  the  third  task.  It 
contains  number  of  loads  connected  to  a  certain  bus  and  max¬ 
imum  demand  of  every  load. 

7.2.  Assessment  of  branch  power  flow  and  branch  loss 

Assessments  of  branch  power  flow  and  branch  loss  are  con¬ 
ducted  in  this  step.  This  should  be  done  in  order  to  know  whether 
or  not  power  flow  on  the  transmission  line  violates  the  maximum 
capacity  of  the  branch.  If  the  power  flow  is  under  the  range  of 
maximum  capacity  then  the  system  will  carry  out  uniform  pricing 
stipulation.  Once  it  exceeds  the  thermal  limits  a  re-dispatch  process 
is  required  to  ensure  that  power  system  remains  secured  for 
transactions  while  keeping  a  balanced  supply  and  demand  function. 

7.3.  Modeling  of  IMP  schemes 

When  redispatch  is  necessary,  different  model  of  LMP  are  for¬ 
mulated.  Three  schemes  of  the  LMP  are  introduced:  LMP-lossless, 
LMP-loss  and  LMP-TUT. 

7.3.1.  Scheme  I:  Named  as  LMP-lossless  model,  has  two  components 
which  are  energy  price  (EP)  and  transmission  congestion  revenue 
(CR),  while  transmission  losses  are  ignored 

LMPiossiess  =  EP+CR  (8) 

Objective  function  and  its  constraints  are  combined  to  formu¬ 
late  a  Lagrange  function  using  Lagrangian  multipliers.  These 
multipliers  are  referred  to  as  dual  variables  or  shadow  prices: 

nc  no  /  \ 

UPDj,  P(,i,  h.  Pi)  =  Z  Ci(Pc,i)  -  Z  Bj(Pn jO+Aj  Z  Pd)  Z  Pci  ) 
i=l  j=l  \j  i  J 

—  Z  Pmin.GilPa  -  Pci " )  +  Z  Pmax.Ci^Ci  ~  PcT ) 
i  i 

-  Z  PminJlow,m(fl°Wm  -fl0W™n ) 
m 

+  Z  Pmaxflow,m(fl°Wm-floW™ax)  (9) 


LMPi0SS  =  EP+CR+LC  (11) 

Formulation  of  network  losses: 

PdPa)  =  \flowm]T[R]\flowm] 

[flowm]  =  Z  SmiPGi  —  Z  SmjPDj  (12) 

i  j 

The  branch  flow  may  be  stated  in  a  general  form: 
flow  =  S.Pj  (13) 

Hence 

Pl  (Pa)  =  [Pi]r[S]'  [R][S][P,]  =  ’[Pi]'  [G,j[P,] 

For  [Gt]  =  2[S]t[R][S]  (14) 

Both  Sm,(i)  and  Sm  y)  are  elements  of  the  shift  factor  related  to  the 
network  connection  node  between  the  generator  or  supplier  i  and 
load  or  customer  j.  The  Lagrange  function  is  then  formulated  as: 

nc  nD 

L(PDj,Pc,iJi,Pi)  =  Z  Ci(Pci)  -  Z  Bj (Pof) 

i  =  t  1=1 

+  Z  PDJ  -  Z  Pa  -  PdPa)  )  -  Z  Pmin,Dj(PDj  ~  Pp f ) 

+  Z  Pmax.Dj(Pp i  -  Pd“X)  -  Z  Pminfii(Pa  ~  Pc! ") 

i  > 

+  Z  Pmax,C,i(Pci  -  Pci*)  ~  Z  PminJlow,m(fl°wm 
i  m 

-flow™" )  +  2  Pmax/low.m(floWm  -flow™* )  (15) 

m 

The  LMPs  may  be  derived  as  the  first  order  KKT  optimality 
condition  as  follows: 

Pi=Aj<I—  Z  Gj.n_n'|  Z  Pd —  Z  Pd,)  1  +Z Pmmjl0W,m  Sm 
l  i'  =  l  \iel(()  jeJti')  )  J  m 

—  Z  Pmaxflow.m  ■~’m  (16) 

m 

7.3.3.  Scheme  III:  Named  as  LMP-TUT  is  formulated  based  on  the 
LMP-loss  but  takes  into  account  a  tariff  for  transmission  usage  (TUT) 
as  well 

LMPtut  =  EP+CR+LC+TUT  (17) 

Total  transmission  usage  tariff  is  written  as: 

Pm(Pa)  =  OCmfloWm  (18) 

Formulation  of  network  losses 


578 


M.  Bachtiar  Nappu  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  38  (2014)  572-580 


The  model  for  network  losses  refers  to  Eq.  (12).  The  Lagrange 
function  becomes: 


+*i  Z  Pdj  -  Z  Pa-  PdPa )  -  Z  Vmin,Dj(PDJ  ~  Pd™  ) 


P(P Dj ■>  P Gi?  5  f^i) 


nG  nD 

Tm(Pa)+  2  Ci(PGi)  ~  2  BjiPDj) 

l  i  —  1  j  =  1 


+  2  AWD/^D/  -  Ptf*)  -  2  Pmin,GiiPci  ~ 
j  i 
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Fig.  3.  The  flowchart  of  LMP  Schemes  Improvement. 
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+  Z  Pmax,Gi(PC.i  ~  Pc/“)  “  Z  Mmin/low,m(f^Wm  -//0W™'n ) 
i  m 

+  'LVmaxJlow,m(fl°Wm  -//OVV^”) 

m 


(19) 


Thus,  the  suppliers’  and  customers’  LMPs  are  stated  as: 


{ =  i 


Pi~  hi  1  -  Z  Gin-nl  Z  PGi  Z  PDj  )  f  +  Z  Pminjlow.m  -S‘ 


i  e  I(i')  j  ej(j') 


—  Yu  S  —  y  a  h0Win(P(,i)  r 


(20) 


where  /?j  is  locational  marginal  price  at  bus  i,  A,p  are  Lagrange 
multipliers,  i  is  supplier  index  (1,2,3,  ...,  /),  j  is  consumer  index 
(1,2,3,  S  is  shift  factor  of  the  transmission  network,  m  and  n 
are  line  number  where  the  supplier  i  or  customer  j  is  connected,  g, 
is  output  of  generator  i,  dj  is  demand  of  load  j,flowm  is  power  flow 
on  line  m,  nG,  nD  are  number  of  generators  and  load,  CL  is 
transmission  network  loss  representation  matrix,  am  is  transmis¬ 
sion  charge  rate  for  line  m  ($/MW  h),  Q  is  production  cost  at  bus  i, 
Bt  is  consumer  benefit  at  bus  i,  Pa  is  active  power  output  of  ith 
generator  and  PDj-  is  active  power  demand  of  consumer  j. 


7.4.  Developing  optimisation  tools 


impacts  such  as  market  inefficiency  owing  to  maximum  capacity 
of  the  network  and  the  possibility  of  exercising  market  power  by 
strategic  behaviour  of  the  producers  which  may  create  high  prices 
in  a  congested  area. 

DC-OPF  has  been  widely  used  for  dispatching  power  and 
clearing  energy  to  determine  the  LMP  due  to  its  speed  and 
robustness,  particularly  in  market  simulation  and  planning.  How¬ 
ever  DC-OPF  does  not  consider  the  network  losses.  As  the  DC-OPF 
is  modified  by  incorporating  the  line  losses,  such  advantageous 
features  do  not  exist  anymore  due  to  the  complexity.  Therefore  a 
new  approach  of  optimal  power  flow  to  maintain  the  linearity  and 
superposition  features  of  the  LMPs  model  is  proposed  while  it  is 
still  able  to  account  for  both  congestion  cost  and  losses  cost.  The 
proposed  method  uses  shift  factor  based  optimal  power  flow. 
Optimisation  scenarios  developed  in  this  research  is  based  on  the 
perfect  competition  where  each  market  participant  is  performed 
by  their  marginal  cost.  The  shift  factor  is  the  actual  component 
that  determines  the  power  flow  over  a  given  transmission  line 
from  the  source  node  (generation)  to  the  sink  node  (load)  which 
can  be  obtained  through  three  stages  based  on  a  DC  load  flow 
network  model. 


After  creating  an  LMP  model  for  three  schemes  optimisation  tools 
are  developed.  With  reference  to  their  specific  models  tools  are  made 
in  such  a  way  that  they  are  able  to  examine  and  evaluate  the  optimal 
condition  of  every  scheme  through  an  iterative  process. 

7.5.  Examining  the  most  efficient  reference  node 

As  shift  factor  is  a  reference  node  based  optimal  power  flow 
methodology,  used  for  a  particular  node  that  can  obtain  the  lowest 
objective  value.  Once  an  optimal  condition  is  reached,  the  objec¬ 
tive  value  of  the  schemes,  generator's  output  and  total  branch 
losses  are  calculated.  Since  market  settlements  are  cleared  at  their 
related  LMPs,  the  market  clearing  price  will  yield  a  residue  which 
is  known  as  merchandising  surplus.  Rudnick  et  al.  [74]  referred  to 
this  term  as  transmission  revenue,  while  [75]  named  it  as  a 
marginal  network  revenue.  By  computing  LMPs  at  all  nodes 
merchandising  surplus  is  analysed  to  get  transmission  congestion 
revenue,  transmission  losses  cost  and  transmission  usage  tariff. 
These  value  based  optimisation  processes  are  beneficial  for  further 
assessment  of  work  plans  such  as:  analysing  market  power  [76], 
optimising  charge  of  transmission  services  [77],  evaluating  co¬ 
optimisation  of  energy  and  spinning  reserve  [78],  investigation  of 
contingency  screening  [79]  and  forecasting  short  term  LMP  values 
[72,80,81],  Problem  identification  representing  the  basic  tasks  is 
illustrated  in  the  flowchart  form  of  the  LMP  schemes  improvement 
in  Fig.  3. 

In  order  for  the  proposed  method  above  to  obtain  appropriate 
reference  node  and  to  perform  the  lowest  overall  cost  the  follow¬ 
ing  iteration  process  should  be  undertaken: 

1.  select  any  arbitrary  node  as  reference  prior  to  running  simulation 

2.  choose  particular  node  which  performs  the  lowest  nodal 

price  and 

3.  re-run  the  simulation  after  fixing  that  node  as  the  reference  node. 

8.  Conclusions 

There  are  three  significant  aspects  of  electric  power  networks 
that  should  be  properly  handled  to  achieve  a  transmission  open 
access  environment,  namely:  transmission  congestion  cost,  trans¬ 
mission  losses  cost  and  transmission  usage  tariff.  The  transmission 
congestion  issue  is  the  central  problem  since  it  may  lead  to  further 
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